addpath('npy-matlab-master\npy-matlab')
A1 = readNPY('map_LD_633nm_30K_NiI2_highres_x.npy')
A1=A1*1e4
A1=rot90(A1)
surf(A1,'edgecolor','none')
set(gca,'xlim',[1 61],'ylim',[1 61])
view(2)
%colormap(winter)
caxis([-1.2,1.2])
set(gca,'box','on','Xcolor',[0 0 0],'Ycolor',[0 0 0])
xticks({})
yticks({})
colorbar('Ticks',[-1:0.5:1],'TickLabels',{})
set(gcf,'Position',[0, 0, 380, 300])